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ABSTRACT 

^ \ We refine a technique to measure the absorption corrected ultraviolet (UV) luminosity of star- 

ve) ■ burst galaxies using rest frame UV quantities alone, and apply it to Lyman-limit [/-dropouts at 

z ~ 3 found in the Hubble Deep field (HDF). The method is based on an observed correlation 
between the ratio of far infrared (FIR) to UV fluxes with spectral slope (3 (a UV color). A 
simple fit to this relation allows the UV flux absorbed by dust and reprocessed to the FIR to be 
| calculated, and hence the dust-free UV luminosity to be determined. International Ultraviolet 

Explorer spectra and InfraRed Astronomical Satellite fluxes of local starbursts are used to cal- 
ibrate the Fpm/ Fiqqo versus (3 relation in terms of Ai6oo (the dust absorption at 1600A), and 

■ the transformation from broad band photometric color to (3. Both calibrations are almost com- 
', pletely independent of theoretical stellar population models. We show that the recent marginal 

""^5 ' and non-detections of HDF [/-dropouts at radio and sub-mm wavelengths are consistent with 

their assumed starburst nature, and our calculated ^4i600- This is also true of recent observa- 
tions of the ratio of optical emission line flux to UV flux density in the brightest [/-dropouts. 
This latter ratio turns out not to be a good indicator of dust extinction. In [/-dropouts, ab- 
solute magnitude Mi6oo,o correlates with (3: brighter galaxies are redder, as is observed to 
be the case for local starburst galaxies. This suggests that a mass-metallicity relationship is 
j> \ already in place at z 3. The absorption-corrected UV luminosity function of [/-dropouts 

■ extends up to Mi6oo,o ~ — 24 ABmag corresponding to a star formation rate ~ 200.M©yr _1 
(Hq = 50 kms" 1 Mpc" 1 , qo = 0.5 assumed throughout). The absorption-corrected UV lumi- 
nosity density at z m 3 is /3i600,0 > 1-4 x 10 ergs -1 Hz -1 Mpc~ 3 . It is still a lower limit since 
completeness corrections have not been done and because only galaxies with A1600 ^3.6 mag 
are blue enough in the UV to be selected as [/-dropouts. The luminosity weighted mean dust 
absorption factor of our sample is 5.4 ± 0.9 at 1600A. 

Subject headings: galaxies: starburst - early universe - ultraviolet: galaxies - infrared: 
galaxies - radio continuum: galaxies 

1. INTRODUCTION ing capabilities, especially in the optical with the 

In the past few years we have seen the fruition of Hu * ble S P ace Tdesc °P e (Williams et al. 1996), (2) 
a great astronomical endevour: the measurement of muto-object spectroscopy, especially down to V « 26 

the rate of cosmological evolution via direct censuses TlQQfi 1 b L 6 "n l^ 1 ^QQtT d (3^ h 
of star formation tracers, metal content, and metal e ' a ' ' owen ^ a e ^ ■• )> an ( ) e 
production rate (Lilly et al. 1995; Pei & Fall 1995; refinement of photometric redshift techniques (e.g. 
Lilly et al. 1996; Madau et al. 1996; Connolly et al. Connolly et al 1997) especially the Lyman-hmit 
1997; Madau, Pozzetti, & Dickinson 1998). This technique (e.g. Steide & Hamilton 1992; Steidel, Pet- 
achievement has come about thanks to several key tim > & Hamilton 1995) which allow redshift estimates 
developments including (1) deep sub-arcsecond imag- 

1 Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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to be made to greater depth, and out to z « 5. 

In a pioneering paper Madau et al. (1996; hereafter 
M96) derived the z evolution of the metal produc- 
tion rate pz(z) (in comoving coordinates), which is 
directly proportional to the volume density of the 
star formation rate (SFR). The "Madau plot" pro- 
vides a direct constraint on the evolution of galaxies, 
hence it's importance. It shows a strong (factor ~ 
10) decline in pz(z) from its peak at z ~ 1 — 2 to the 
present epoch, as pointed out earlier by Lilly et al. 
(1996). This decline mirrors that seen in the lumi- 
nosity density of quasars (Boyle & Terlevich 1997). 
On the other side of the peak, at 2 » 4, Pz{z) is 
apparently only equal to that of the present epoch. 
Subsequent work by Connolly et al. (1997) fills in the 
z = 1 — 2 gap and is consistent with M96. The over- 
all shape of the curve is consistent with models of 
the observed evolution in the properties of damped 
Lyman-a systems (Pei &: Fall, 1995). This form for 
pz( z ) is consistent with hierarchical models of galaxy 
formation (e.g. Searle & Zinn, 1978; Zepf & Ashman, 
1993; Baugh et al. 1998). 

One concern about the Madau plot is that it is 
heavily biased to observations in the rest frame UV. 
The problem is that dust efficiently scatters and ab- 
sorbs UV radiation, greatly complicating the inter- 
pretation of the detected UV emission. Not only 
is the amount of dust important, but so is the ge- 
ometry (where it is relative to the stars), and the 
composition of the dust. For example, in the Local 
Group the UV "extinction law" varies from galaxy to 
galaxy and within galaxies, (e.g. Fitzpatrick 1986). 
The term "extinction law" was originally defined for 
stars. It quantifies the wavelength dependent com- 
bined effect of dust absorption and scattering out of 
the line of sight toward point sources. The geometry 
is quite different in distant galaxies where individ- 
ual stars can not be isolated, and instead galaxies 
in their entirety, or large portions of them are ob- 
served. In such cases much of the scattered light may 
be recovered (by scattering into the line of sight), 
and it matters a great deal how porous the dust dis- 
tribution is and whether the stars are mixed with the 
dust or not (e.g. Witt, Thronson, & Capuano 1992; 
Calzetti, Kinney & Storchi-Bergmann 1994). There 
are a variety of terms for the wavelength dependent 
net dimming of the light including "effective extinc- 
tion law" (Calzetti et al. 1994) and "obscuration law" 
(Calzetti 1997). Here we adopt the term "absorption 
law" to emphasize that it is dust absorbing photons 
that causes the dimming. Absorption laws can be 
derived empirically (e.g. Calzetti et al. 1994; Calzetti 
1997) and hence average together all the above ge- 
ometric factors into a simple wavelength dependent 
correction. 



Observations of high redshift galaxies in the sub- 
mm and mm bands prove that molecular gas and 
dust formation have occurred as far back as z = 4.7 
(Omont et al. 1996; Ivison et al. 1998; Guilloteau 
et al. 1997; Frayer et al. 1999). Especially interesting 
are the recent spate of blank field sub-mm detection 



with SCUBA, which we discuss in more detail in §5.2 
and §^3. These can account for a significant fraction 
of the sub-mm background, and are mostly attributed 
to dusty star forming galaxies having z > 1. Follow- 
ing from UV observations of local starbursts (Meurer 
et al. 1995, hereafter Paper I) we have argued that 
the correction factor to pz{z) is considerable, a factor 
~ 15 at z ~ 3 (Meurer et al., 1997; hereafter Paper 
II), while others (e.g. Madau 1997; Pettini et al. 1997) 
prefer fairly modest, ~ factor of three, corrections. 
The amount of high-z absorption has a direct bearing 
on the interpretation of the Madau plot; large correc- 
tions at high-z can push back the peak pz, as might 
be expected by monolithic collapse dominated galaxy 
formation (e.g. Eggen, Lynden-Bell, & Sandage 1962; 
Bower, Lucey & Ellis 1992; Ortolani et al. 1995). 

Here we reexamine the UV luminosity density at 
high-z, with particular emphasis on the [/-dropouts 
in the HDF. Our technique has two strengths. Firstly, 
it is based on the similarity between local starburst 
galaxies and Lyman limit systems as summarized 
in § ^| The main assumption we make is that the 
high redshift galaxies have the same multi-wavelength 
spectral properties as local starbursts. We use lo- 
cal starburst galaxies to calibrate the relationship 
between UV reddening and UV absorption (in § ||) 
thus largely avoiding the uncertainties in population 
model fitting and the choice of effective UV absorp- 
tion laws. Secondly, our calibration of the UV ab- 
sorption law recovers the intrinsic (i.e. dust-free) UV 
emission that is absorbed by dust and reradiated in 
the far infrared. Hence, up to the modest amounts of 
absorption considered here (a few magnitudes), UV 
observations are capable of recovering the total bolo- 
metric luminosity of starbursts. We apply our cali- 
brations to the [/-dropouts in the HDF in § ||. In § [| 
we demonstrate the validity of our starburst assump- 
tion by showing that new constraints on the SEDs 
of these high-z systems are consistent with the SEDs 
of local UV selected starbursts. Implications of our 
results are discussed in §||, while § summarizes the 
paper. 

Throughout this paper we adopt Hq = 
SOkms^Mpc" 1 , q = 0.5. 



2. HIGH-Z GALAXIES AND LOCAL STARBURSTS 

Progress in understanding the Lyman-dropout 
galaxies would be accelerated if we could find analogs 
to these galaxies in the local universe. Fortunately, 
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there is mounting evidence that the Lyman dropout 
galaxies indeed resemble starburst galaxies in many 
respects. 

While the magnitude of the extinction corrections 
to the rest-frame UV light (the subject of the present 
paper) is still a matter of debate, it is nevertheless 
clear that local starbursts have surface-brightnesses 
that are similar to the distant galaxies (Paper II). 
That is, the high-z galaxies appear to be 'scaled-up' 
versions of the local starbursts with the same SFRs 
per unit area and comparable stellar surface-mass- 
densities. In this regard, starbursts are a far better 
match to the Lyman-dropout galaxies than are or- 
dinary late-type galaxies, which have SFRs per unit 
area that are several orders-of-magnitude smaller (cf. 
Paper II; Kennicutt 1998). 

The overall rest-frame UV-through-optical spectral 
properties of local starbursts and the Lyman dropout 
galaxies are similar in terms of the spectral energy 
distribution (cf. Sawicki & Yee 1998) and the dom- 
inance of the UV and optical spectra by absorption 
and emission lines respectively (cf. Leitherer 1997; 
Pettini et al. 1998). The Lyman limit systems in the 
HDF are clearly spatially resolved (Giavalisco, Stei- 
del, & Macchetto 1996), demonstrating that their rest 
frame UV light is not dominated by an AGN. In ad- 
dition the spectra of the Lyman-dropout galaxies is 
more similar to starburst galaxies than AGN. Indeed 
the typical strengths of the UV stellar and interstellar 
absorption-line features are comparable in starbursts 
and Lyman-dropout galaxies, indicative of broadly 
similar stellar and interstellar content (e.g. Conti, 
Leitherer, & Vacca 1996; Lowenthal et al. 1997). 
The local starbursts and Lyman dropout galaxies also 
have similar distributions of observed UV colors, (Pa- 
per II, and here in §Q), - which we argue chiefly mea- 
sures the amount of dust-reddening. 

Even the dynamical state of the interstellar 
medium in the Lyman dropout galaxies resembles 
that in local starbursts. Both are driving extensive 
outflows of gas at velocities of several hundred km 
s -1 . The signature of this is the systematic blueshift 
of the broad UV interstellar absorption-lines with re- 
spect to the rest-frame optical (e.g. H/3) and/or rest- 
frame UV (Lya) nebular emission-lines (Pettini et al. 
1998; Franx et al. 1997; Kunth et al. 1998; Gonzalez- 
Delgado et al. 1998). The velocity dispersion of the 
emission-line gas is also similar in the two galaxy 
classes (Pettini et al. 1998). 

To reiterate, in all their properties measured to- 
date, the Lyman dropout galaxies appear to be sim- 
ply scaled-up versions of local starbursts: larger and 
more luminous, but otherwise remarkably similar. 
The challenge is to determine the intrinsic luminosi- 
ties from the UV fluxes and colors of dusty Lyman 



limit systems. The above similarities justifies our ap- 
proach of calibrating the effects of dust on local star- 
bursts and applying it to the high-z galaxies. 

3. THE LOCAL SAMPLE CALIBRATION 

Figure [l] illustrates the heart of the technique. It 
shows the relationship between the ratio of far in- 
frared (FIR) and ultraviolet (UV) fluxes (here at 
1600A) and the UV spectral slope (3 (equivalent to 
a color) for a sample of starburst galaxies. Since the 
FIR flux is due to dust radiatively heated by the ab- 
sorbed UV radiation, the y axis is a measurement of 
dust absorption. This figure then shows that dust 
absorption is correlated with UV reddening. Such a 
relationship is expected for dust predominantly lo- 
cated in a foreground screen (Witt et al. 1992; Paper 
I), although this screen need not be homogeneous 
(Calzetti et al. 1994; Calzetti 1997). Because this 
correlation links -Ffir to UV quantities, then regard- 
less of the exact dust geometry, this figure provides 
a powerful empirical tool for recovering the radiation 
reprocessed by dust and thus determining the total 
absorption-corrected UV flux of starbursts, using UV 
quantities alone. 

In the next sub-section we define our local cali- 
brating samples and the quatities we use. In order 
to apply this tool and interpret the results we per- 
form three calibrations. 1. -Ffir/^igoo is calibrated 
in terms of absorption at 1600A and this is fitted 
clS cl function of (3. 2. The spectroscopic index j3 
is calibrated in terms of broad-band colors (with a 
z correction), since photometric colors are easier to 
measure than spectroscopic ones. 3. The relationship 
between luminosity measured at 1600A is related to 
SFR. These calibrations are detailed in the final three 
sub-sections. 



3.1. Sample and definitions 

The local sample used to derive the various calibra- 
tions in this paper is listed in Table [I]. It includes the 
sources shown in Fig. EL and is drawn from the Inter- 
national Ultraviolet Explorer (IUE) atlas of Kinney 
et al. (1993), and contains galaxies having an "ac- 
tivity class" consistent with being a starburst, i.e. 
SB nuc. (starburst nucleus), SB ring (starburst ring), 
BCDG (blue compact dwarf galaxy), or BCG (blue 
compact galaxy). Although the UV sources in these 
galaxies tend to be centrally concentrated we find 
that galaxies with optical diameters -D25 > 4' (D25 
is measured at B = 25 mag arcsec -2 ) tend to fall 
above the relationship shown in Fig. |]. This is prob- 
ably because significant UV emission extends beyond 
the 20" x 10" IUE aperture. Hence these large galax- 
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Fig. 1. — The ratio of FIR to UV flux at I600A compared to UV spectral slope f3 for UV selected starburst galaxies. The right 
axis converts the flux ratio to 1600A absorption Aieoo using eq. M. The solid line shows our adopted linear fit to the ^1600, (3 
relationship. The dotted line shows the proposed dust absorption/population model of Pettini et al. (1997). 



ies were excluded from the sample 2 . Further limiting 
the sample to galaxies having D25 < 2.5' removes a 
few more outliers, but not just galaxies above the re- 
lationship. Furthermore, it also severely depletes the 
points with (3 > —0.5. Applying fits to data limited 
in this way changes our final pimo(z = 2.75) results 
by < 10%, hence we retain D25 < 4' as the diame- 
ter limit for the local sample. The lack of systematic 
residuals for galaxies up to 12 times larger than the 
IUE aperture indicates that the UV emission of these 
galaxies is very compact, e.g. in a circum-nuclear star- 
burst. 

The IUE spectra in the Kinney et al. atlas were 
measured to determine the UV quantities required 
for Fig. [j]. The ultraviolet flux at 1600A, F1600 is a 
generalized flux of the form F\ = \f\ , and fx is the 
flux density per wavelength interval. It was measured 
with the IRAF / STSDAS 3 package SYNPHOT em- 
ploying a square passband with a central wavelength 
of 1600A and width of 350A. This filter is meant to 
approximate the rest frame parameters of the stan- 
dard WFPC2 filters F606W and F814W for objects 



with redshifts z = 2.75, and 4, respectively (i.e. U and 
B band dropouts). The ultraviolet spectral slope f3 
is determined from a power law fit of the form 



fx oc X 13 . 



(1) 



to the UV continuum spectrum as defined by the ten 
(rest wavelength) continuum bands listed by Calzetti 
et al. (1994). These spectral fits were performed af- 
ter first removing Galactic extinction using the law 
of Seaton (1979) and taking Galactic extinction val- 
ues Ab = A.\E{B — V) from Burstein & Heiles 
(1982, 1984) as listed by NED 4 . Since the contin- 
uum spectrum is never exactly a pure power law, 
(3 is subject to systematic uncertainties due to the 
placement of the continuum windows. Eighteen of 
the data points in Fig |l| represent galaxies observed 
with one of IUE's short wavelength (SW) cameras 
(A ~ llOOA to 1975A) only and not with either of the 
LW cameras (A « 1975A to 3000A), hence they do 
not have data in Calzetti's et al. (1994) tenth window 
(A = 2400 - 2580A). Their j3 values were determined 
from the SW only measurements using the following 

2 The other galaxies that were excluded are NGC 1569 because of excessive foreground Galactic extinction; NGC 3690 because 
the IUE pointing is likely to be wrong (Paper I); the BCDGs Mrk209, Mrk220, and Mrk499 because they have neither IRAS fluxes 
nor upper limits; and finally five galaxies with low S/N IUE spectra in the Kinney et al. atlas: NGC 4853, IC 2184, Mrk 309, 
Mrk 789, UGC 6448. 

3 The Image Reduction and Analysis Facility software package is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. STSDAS is the Space Telescope Science Data Analysis Software package for IRAF, distributed by 
STScI. 

4 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space Administration. 
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Table 1 
Local sample: relevant data 



name 


A r (Gall 


8 


81 F 1600 ' 




1ok( Fh/3 ) 

Sl /l600 > 




F [OIII]5007 
Fk/3 


Relevant 




[mag] 






[erg/cm 2 /s] 


[A] 






figures 


NGC 4861 


0.00 


-2.46 


-0.08 


-9.97 








3 


I Zw 18 


0.02 


-2.43 


< -0.84 


< -11.50 








1,3 


NGC 1705 


0.18 


-2.42 


-0.76 


-10.29 


-0.01 


2.85 


4.61 


1,7 


Mrk 153 


0.00 


-2.41 


-0.55 


— 10.92 








1 


Toll924-416 


0.30 


—2.12 


—0.12 


— 10.17 


0.94 


2.91 


4.87 


1,3,7 


UGC 9560 


0.00 


-2.02 


-0.03 


— 10.41 


0.79 


3.36 


4.13 


1,3,7 


Mrk 66 


0.01 


— 1.94 


0.25 


— 10.56 


0.68 


2.85 


2.90 


1,7 


NGC 3991 


0.02 


—1.91 


0.34 


-9.80 








1,3 


NGC 3738 


0.00 


-1.89 












3 


UGCA 410 


0.02 


-1.84 


-0.07 


-10.86 








1 


Mrk 357 


0.17 


-1.80 


0.21 


— 10.37 


0.80 


3.25 




1,7 


NGC 3353 


0.00 


-1.79 


0.65 


-9.61 








1,3 


MRK 54 


0.02 


-1.78 


0.28 


-10.27 








1 


NGC 1140 


0.11 


-1.78 


0.32 


-9.77 


0.96 


3.18 


2.66 


1,7 


Mrk 36 


0.00 


— 1.72 


—0.26 


— 10.94 








1 


MCG6-28-44 


0.02 


-1.77 


0.31 


-10.27 








1 


NGC 1510 


0.00 


—1.71 


0.10 


-10.36 


0.67 


3.11 


4.95 


1,3,7 


MRK 19 


0.14 


-1.69 


0.10 


-10.64 








1 


NGC 4214 


0.00 


-1.69 












3 


NGC 4670 


0.04 


-1.65 


0.17 


-9.85 








1,3 


NGC 1800 


0.00 


-1.65 






0.41 


3.08 


3.99 


7 


UGC 5720 


0.00 


— 1.62 


0.22 


-9.97 








1,3 


UGC 5408 


0.00 


-1.60 


0.56 


— 10.22 








1 


NGC 7673 


0.14 


-1.50 






0.75 


4.46 


2.13 


3,7 


NGC 3125 


0.25 


-1.49 


0.55 


-9.64 


1.07 


3.28 


4.72 


1,7 


Haro 15 


0.08 


-1.48 


0.31 


— 10.16 


0.74 


2.85 


2.32 


1,7 


NGC2537 


0.15 


— 1.44 


1.05 


-9.73 








1 


UGC 3838 


0.07 


—1.41 


0.26 


-10.55 








1 


NGC 7793 


0.00 


— 1.34 






0.63 


4.60 


0.55 


3,7 


NGC 5253 


0.18 


-1.33 


0.69 


-8.86 


1.02 


3.18 


4.68 


3,7 


NGC 7250 


0.62 


-1.33 


0.46 


-9.77 


0.80 


3.14 


3.48 


1,7 


Mrk 542 


0.07 


— 1.32 






0.54 


4.00 


0.79 


7 


NGC 7714 


0.15 


— 1.23 


0.84 


-9.32 


1.18 


4.36 


1.60 


1,3,7 


Mrk 487 


0.02 


—1.15 






1.01 


3.25 


6.61 


7 


NGC 3049 


0.04 


—1.14 


0.85 


-9.84 


0.99 


4.00 


0.23 


1,7 


UGC 6456 


0.09 


—1.10 


0.11 


— 10.74 








1 


NGC 3310 


0.00 


-1.05 


1.01 


-8.83 








1 


NGC 5996 


0.06 


-1.04 


0.98 


-9.64 


0.86 


4.76 


0.38 


1,7 


NGC 4385 


0.03 


-1.02 


0.97 


-9.65 


0.97 


5.48 


0.79 


1,7 


Mrk 499 


0.00 


-1.02 






0.38 


4.65 


3.57 


7 


NGC 5860 


0.03 


-0.91 


0.92 


-10.04 


0.88 


6.04 


0.64 


1,7 


IC 1586 


0.07 


-0.91 


0.71 


-10.28 


0.89 


5.30 


1.35 


7 


NGC 2782 


0.00 


-0.90 


1.17 


-9.33 








1 


ES0383-44 


0.13 


-0.83 


0.71 


-10.13 








1 


NGC 5236 


0.14 


-0.83 






0.69 


3.91 


0.22 


7 


NGC 2415 


0.16 


-0.80 


0.97 


-9.36 








1 


MCG-01-30-33 


0.05 


-0.77 


1.09 


-9.84 








1 


NGC 1614 


0.22 


-0.76 


2.07 


-8.84 


1.24 


7.85 


0.96 


7 


NGC 6217 


0.14 


-0.74 


1.26 


-9.22 


0.79 


5.08 


0.39 


1,7 


NGC 6052 


0.10 


-0.72 


1.14 


-9.48 


1.04 


3.54 


2.01 


1,7 


NGC 4500 


0.00 


-0.67 


1.13 


-9.69 








1 


IC 214 


0.11 


-0.61 


1.35 


-9.56 


0.66 


5.08 


1.82 


7 


NGC 3504 


0.02 


-0.56 


1.45 


-8.96 








1 


NGC 4194 


0.00 


-0.26 


1.63 


-8.99 


1.24 


6.81 


1.13 


1,7 


NGC 2798 


0.01 


0.03 


2.05 


-9.01 








1 


NGC 3256 


0.55 


0.16 


1.88 


-8.36 


1.23 


5.48 


0.47 


1,7 


NGC 7552 


0.00 


0.48 


2.12 


-8.44 


1.29 


6.11 


0.1 


1,3,7 
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relationship: 



{P ~ Psw) 



-0.16 ±0.04. 



(2) 



This was determined from measuring (3 of 16 
high signal-to-noise IUE spectra with no noticeable 
SW/LW break, both with and without the tenth win- 
dow. The uncertainty is the standard error on the 
mean. 

The only non-UV quantity in Fig. [l] is the far- 
infrared flux -Ffir which is derived from InfraRed As- 
tronomical Satellite (IRAS) 60 and 100/mi flux den- 
sities listed by NED and is defined by Helou et al. 
(1988). For convenience in notation we define the 
ratio F F i R /Fi 600 = IRXi 600 

3.2. Errors 

Typical error bars are shown in the bottom right of 
fig |l[ We estimate an uncertainty in (3 of 0.14 from 
the x 2 fits t° the spectra, and the scatter about the 
mean relationship given in eq. ||). We adopt an av- 
erage uncertainty in Fi6oo of 5%, which we showed 
in Paper I is the typical level of agreement between 
broad band HST UV fluxes and IUE spectra. This is 
higher than expected from the error spectra presented 
with the Kinney et al. (1993) IUE atlas, which are 
commonly at the 10% - 20% level per pixel, and hence 
< 1% when integrated over the UV passband. The 
level of agreement between HST and IUE is worse 
probably because of the centering errors with IUE. 
The median uncertainty in Ffir for our sample is 
7%, derived from the uncertainties in Fqq and F100 as 
listed by NED. Combining the -F1600 and -Ffir errors 
we estimate that the error in log(FFiR/Fi6oo) is 0.04 
dex. 

These uncertainties represent random errors. Sys- 
tematic errors will effectively move all data points 
with respect to the zeropoints of the axes. C. Robert 
(1999, private communication) estimates that there is 
a sytematic uncertainty of 0.15 in (3 due to the place- 
ment of the continuum windows in the spectrum fit. 

Figure |l| illustrates that the random errors in the 
quantities are smaller than the scatter about the rela- 
tionship. Hence, the scatter in Fig. |] is largely intrin- 
sic. Possible sources of intrinsic scatter are dust ge- 
ometry, dust composition, and intrinsic spectral slope 
(which depends on the star formation history), all of 
which may vary from galaxy to galaxy. 

3.3. (3 - Aiqqq calibration 

We assume that -Ffir is due to the thermal emission 
of dust heated by the radiation it absorbs. This con- 
sists primarily of non-ionizing photons with the ad- 
dition of Lya photons which are resonantly scattered 

5 Starburst population models indicate that typically only < 
912A (Leitherer & Heckman, 1995). The fraction of this that is d 
Biermann, & Mezger 1978). 



by Hydrogen until they are absorbed by dust grains. 
We assume that ionizing photons do not significantly 
contribute to dust heating 5 . For these assumptions 
IRX1600 can be written as 



IRXi 



-F1600 



%. + J 9 7 2X /v,o(l-10- 



fl600,0l0-°- 4A 1600 



F b oi 



(3) 



fx,o is the unattenuated flux density of the emitted 
spectrum, A\ is the net absorption in magnitudes by 
dust function of wavelength, and FL yQ , is the 

Lya flux. In our models FL yQ , is derived from the 
spectrum shortward of 912A using the standard as- 
sumption that each ionizing photon results in one re- 
combination and consequent decay down to case-B 
population levels. Customarily one adopts a form of 
the absorption (or extinction) law scaled by the op- 
tical reddening, e.g. A\ = X\E(B — V). In eq. |3|the 
first term gives the bolometric flux of the absorbed 
radiation, normalized by F1600! and the second term 
gives the fraction of the bolometric flux emitted by 
dust that is intercepted by the FIR passband (e.g. 
Helou et al. 1988). 

Strictly speaking IRX1600 i n our model depends 
on the intrinsic spectrum and the absorption curve, 
even at wavelengths well outside the 1600A passband. 
However, for dust heating dominated by young stellar 
populations eq. |3]is well approximated by the formula 



IRXieoo = (l0°' 4Al e°° - l) 



-Fl600,0 



V F bol J 



Dust 
(4) 

Here the first term gives the fraction of Fi6oo,o that 
is absorbed by dust; the second term is similar to a 
bolometric correction, giving the maximum amount 
of heating available to the dust divided by the in- 
trinsic 1600A flux; and the third term is the same as 
the second term in eq. ^. The ratio of eqs. ||| and || 
asymptotes to unity as A\ — > 00 for all A. We have 
calculated this ratio under a variety of assumptions 
applicable to starbursts. The constant star formation 
rate (CSFR) solar metallicity models of Leitherer and 
Heckman (1995) having a Salpeter (1955) IMF slope 
and durations of At = 1 to 100 Myr were used as the 
initial spectra. Leitherer and Heckman (1995) varied 
other parameters in their population models includ- 
ing metallicity, IMF slope and cutoff, and other star 
formation histories (i.e. instantaneous bursts). They 
find that as long as the population is ionizing (i.e. 
produces a strong ionizing flux, as expected for star- 
bursts by their strong emission line spectra), these 

10% of the bolometric luminosity of starbursts is emitted below 
irectly absorbed by dust grains is estimated to be ~ 25% (Smith, 
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parameters have very little effect on the overall UV 
spectrum which has (3q = —2.0 to —2.6 in all cases. 
This range in intrinsic (3 is well covered by the models 
we examine here. The absorption was parameterized 
by E(B — V) = 0.01 to 1 and a variety of forms of 
X\ including the following extinction laws (assuming 
extinction = absorption): Galactic (Seaton, 1979), 
LMC (Howarth, 1983), and SMC (using a polynomial 
fit to the Bouchet et al. 1985 curve); and the starburst 
absorption laws of Kinney et al. (1994), Calzetti et al. 
(1994), and Calzetti et al. (1997). In all cases the ra- 
tio of equations ||| and |] is within 30% of unity for 
^4i600 ^ 0-3 mag. We conclude that eq. ^ is a good 
approximation of IRX1600 when the dust is heated by 
intrinsically young stellar populations and as long as 
the form of X\ is not extremely different from stan- 
dard absorption/extinction laws. 
Equation can be rewritten as 



IRX 



1600 



(10' 



,0.4Ai 6 oo 



(5) 



where B is the ratio of the two bolometric like cor- 
rections: 

BC(1600)* 

(6) 



B = 



BC(FIR) Dust ' 



BC(1600), 



'F L 



VQ 



+ r i fx od\" 

J 912A _ ,u 



F 



1600,0 



BC(FIR) dust = (P*) 

\£ FIR /Dust 



(7) 



(8) 



Here we assume that B is constant, and estimate the 
best value for it. We do this by cross comparing ob- 
servations with models of the photometric properties 
of stellar populations (BC(1600)*) and dust emission 
(BC(FIR) dust ). 

Values of BC(1600)* were calculated for the CSFR 
population models of Leitherer & Heckman (1995) 
with standard Salpeter (1955) IMF slope, and an up- 
per mass limit of 100 For burst durations in the 
range of to 100 Myr, we find 1.56 < BC(1600)* < 
1.76. We also calculated 60(1600)* directly from the 
observed SED of NGC 1705, one of the least extincted 
sources in our local sample. We started with the 
observed IUE + optical SED published by Storchi- 
Bergmann et al. (1995) after correction for Galactic 
extinction. This was extended down to A = 912A 
by extrapolation using the power law fit (yielding 
(3 = —2.40). The SED was extended redwards to 
A = 2.4 /im using the JHK photometry of Jones as 
reported in Meurer et al. (1992) 6 . The total Ha 
flux given by Meurer et al. (1992) was converted to 
Fhya assuming case B recombination. The integral 
in eq. |^ was evaluated directly from the SED, and 

6 Their I band photometry was used to calculate a 0.34 mag 
Bergmann et al. (1995). 



-^1600.0 was measured with SYNPHOT. This yields 
BC(1600)* = 1.81. Note that this may be a bit 
overestimated (by up to 20%) since there appears to 
be a break between the SW and LW portions of the 
IUE spectrum, with the SW portion being relatively 
underestimated. We conclude that within the un- 
certainties the NGC 1705 measurement is consistent 
with the model estimates. 

The bolometric dust correction BC(FIR)D U st 
was estimated empirically by directly calculating 
BC(FIR)Dust for the ultra-luminous infrared galax- 
ies observed in the FIR-submm by Rigopoulou, 
Lawrence, & Rowan- Robinson (1996). Two galax- 
ies (Mrk 231 and IRAS05189-2524) were rejected 
from their sample because they have strong AGN 
emission. A smooth SED was drawn through the 
SED observations of the remaining sources and in- 
tegrated over 1/xm < A < 1.0mm to determine 
Fboi- The individual BC(FIR)D U st measurements are 
shown in the top panel of Fig. 0, function of 

FIR color / !y (60/im)// I/ (100/xm). For comparison the 
FIR colors of local UV selected starbursts are shown 
in the bottom panel of Fig. [2|. We find that six 
of the seven remaining galaxies in the Rigopoulou 
et al. sample have {BC{FIR)£, us t) = 1.47 with a 
dispersion of 0.13. The seventh, 08572+3915 has 
/ I/ (60/um)//^(100/im) = 1.64, much warmer than the 
majority of the UV selected starbursts, almost all of 
which have 0.3 < / iy (60/im)// iy (100/im) < 1.1. 

The top panel of Fig. |2| shows that the empiri- 
cal BC(FIR)Dust measurements are in good agree- 
ment with Helou's et al. (1988) dust emission models. 
These models are very simple. They treat the dust 
distribution as emitting thermal radiation at a single 
temperature, with emissivity oc v n . We caution that 
more realistic modeling of interstellar dust emission 
includes multiple temperature emission with n ~ 1 
to 2 (e.g. Desert, Boulanger, & Puget 1990). Fig- 
ure H illustrates that single temperature models with 
n ~ 1 provide a good approximation to the dust emis- 
sion of starbursts for the limited purpose of calculat- 
ing BC(FIR)Dust- Adopting this model we have cal- 
culated BC(FIR)Dust for the local UV selected star- 
bursts, which tend to have somewhat cooler colors in 
the FIR than ultraluminous infrared galaxies. These 
have a median BC(FIR)D U st = 1-37, with 80% of the 
sample having 1.36 < BC(FIR) Dust < 1.44. 

We adopt BC(1600)* = 1.66 ± 0.15, and 
BC(FIR) Dust = 1.4 ± 0.2, yielding B = 1.19 ± 0.20. 
Here the "errors" actually indicate the amount of 
leeway allowed by our model calculations and lim- 
ited measurements of starbursts, as described above. 



aperture correction to match the IUE size aperture of Storchi- 
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Fig. 2. — The top panel shows the dust bolometric correction BC(FIR)d ua t for the non-AGN ultra luminous infrared galaxies in 
the sample of Rigopoulou et al. (1996) plotted against FIR color. The dotted, solid, and dashed lines show the single temperature 
dust emission models of Helou et al. (1988), for emissivity indices n — 0, 1, and 2 respectively. The bottom panel shows the 
distribution of FIR colors for local UV selected starbursts. 



Equation || can now be written logarithmically as 

log(IRX 160 o) = log(10 a4Aieo ° - 1) + log B, (9) 

log(IRXi 600 ) = log(10°- 4Al60 ° - 1) + 0.076 ± 0.044. 

(10) 

Here we see the beauty of the approximation eq. it 
allows the direct transformation of IRX1600 to 1600A 
absorption as shown on the right hand side of Fig. |l]. 

Our adopted model for the IRX1600 - P relationship 
is a simple linear fit of vli600 versus (3. To first or- 
der this is what one would expect for absorption by 
a foreground screen of dust. The adopted fit, which 
is shown in Fig. [j], is 

^600 = 4.43 + 1.99^. (11) 

It was determined from separate unweighted least 
squares fits of j4i600 versus 0, and versus ^4i600> 
which were then averaged. The dispersion about this 
fit is 0.55 mag in j4i6oo> 0.28 in and 0.30 dex in 
log(IRXi6oo) (the latter including only galaxies with 
> —2.23). Consequently the standard error in the 
fit zeropoint is 0.08 mag in terms of A1600 an d 0.04 
in terms of 0. Since the scatter about the fit is signif- 
icantly larger than the errors, it is largely intrinsic. 
Hence, a more sophisticated fit (e.g. a full \ 2 tit with 
weighting by errors) will not improve the quality of 
the fit. Instead, improvement can only occur with an 
increase in sample size. 

This fit implies the spectral slope in the absence of 
dust absorption is 0q = —2.23. This result is consis- 



tent with what is expected of naked ionizing popu- 
lations, which Leitherer & Heckman (1995) show to 
have O = -2.0 to -2.6 (cf. Paper I). In Table | the 
slope of the absorption law, dAiQoo/d0 derived here is 
compared with other extinction and absorption laws. 
The slope we derive is in the middle of the range for 
the starburst absorption laws derived in other stud- 
ies. In particular, it is between the values derived by 
Calzetti et al.(1994) and Calzetti (1997). Note that 
the Galactic extinction law has a very large slope be- 
cause the 2 175 A absorption feature counteracts the 
reddening of (3 to a large degree. 

3.4. Photometric versus spectroscopic (3 

In Paper II we derived the following relationship 
between (3 and (Vf306 — -^814)ab color: 

/3phot = 3.23(V(306 — ^814)ab — 2, (12) 

where the color is in the "AB" system (rriAB = 
—48.6 — 2.51og(/^), where f u is the spectral 
flux density per frequency interval in units of 
erg cm -2 s _1 Hz -1 ). This was derived from first prin- 
ciples using the central wavelengths of the F606W 
and F814W WFPC2 filters and assuming a feature- 
less power law spectrum and negligible filter widths. 
However, the UV spectra of starbursts typically con- 
tain numerous absorption features of both stellar and 
interstellar origin (Kinney et al. 1993). Emission 
lines, even Lya, are typically weak in starburst galax- 
ies (Kinney et al., 1993). While the spectroscopically 
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Table 2 
1600a absorption law slopes 



law 


type 


dA 1600 /df3 


reference 


Starburst 


absorption 


1.99 


This work 


Starburst 


absorption 


1.44 


Kinney et al. (1994) 


Starburst 


absorption 


1.97 


Calzetti et al. (1994) 


Starburst 


absorption 


2.30 


Calzetti (1997) 


Galactic 


extinction 


4.42 


Seaton (1979) 


LMC/30 Dor 


extinction 


1.45 


Howarth (1983) 


SMC 


extinction 


1.05 


Bouchet et al. (1985) 



defined (3 avoids the lines, broad band niters cannot. 
Since the concentration of lines is particularly strong 
for 1216 < A < 1700A one may expect an artificial 
reddening of the photometric f3 at z ~ 3 when these 
features are redshifted into the middle of the F606W 
passband. In addition at z > 2.8, Lya enters the 
F606W passband. The Lyman forest strongly sup- 
presses the flux bluewards of Lya in such high redshift 
objects, again reddening /3 pn ot- 

We have used IUE starburst spectra to quantify 
these effects and model the calibration between pho- 
tometric and spectroscopic (3. The spectra we use 
were selected (by eye) to have high S/N and no dis- 
cernible break between LW and SW portion of spec- 
trum. This resulted in a final sample of 15 galaxies, as 
indicated in Table [j]; most of these were also used to 
calibrate the ^4i600 versus (3 fit. Their spectroscopic 
(3 were measured as described above. These spectra 
were redshifted over the range 2 < z < 4. In order to 
adequately account for the F606W filter width, the 
spectrum blueward of Lya was taken to be a pure 
power law (using the measured (3) normalized to the 
mean flux density at A = 1483 ±50 A, and attenuated 
by the Lyman series forest which was modeled follow- 
ing the algorithm of Madau (1995). The (Vesoe — hu) 
colors were then measured using SYNPHOT, yield- 
ing /3 p hot- Figure ^ plots the difference between /3 p hot 
and the standard spectroscopic (3 as a function of z. 

Our calibration of the /3 p hot — (3 relationship is 
a quadratic polynomial fitted to the median points 
shown in Fig. 0: 



which is now on the same scale as the (3 — A\§qq cal- 



ibration of eq. 11 



P phot -13 = 3.22 - 2.66z + 0.545z 2 



(13) 



The fit is a x 2 n t where the errors were taken to be 
the mean difference between the median and the 16th 
and 84th percentiles for a given z. The fit has an rms 
dispersion of 0.06. The reduced x 2 is 0.19, hence the 
model is adequate for the given amount of scatter. 
Combining this with eq. [l^ yields: 

(3 = 3.23(^06 - IbmJab - 5-22 + 2Mz - 0.545z 2 (14) 



3.5. Absolute star formation rate 

The transformation of UV luminosity to SFR is 
highly model dependent. Hence, for most of this pa- 
per the emphasis will be on luminosity rather than 
the SFR. However, for comparison to other studies 
and to get a feel for the physics behind the UV lumi- 
nosity, it is often convenient to transform it to a SFR. 
Our transformation is based on the population mod- 
els of Leitherer & Heckman (1995). We adopt a solar 
metallicity, continuous star formation history popula- 
tion of duration 100 Myr, and a Salpeter (1955) initial 
mass function (IMF) with mass limits of 0.1 to 100 
Mq. For these parameters a SFR of 1 M.q yr _1 pro- 
duces a 1600A flux density of 7.86 x 10 27 erg s^ 1 Hz^ 1 
(9.19 x 10 39 erg s^ 1 A^ 1 ), and hence an absolute mag- 
nitude of -18.15 ABmag (-20.81 STmag); the equiva- 
lent bolometric flux corresponds to 6.4 x 10 9 L@. For 
the same star formation history and IMF parameters, 
Bruzual & Chariot (1993) stellar population models 
produce 8% more flux at 1600A (see also comparison 
in Paper I). Not only is the choice of the particular 
population model important in this calibration, so to 
are the parameters of the IMF and the adopted star 
formation history. For example, if the lower mass 
limit of the IMF is raised to 1 Mq then the luminos- 
ity increases by a factor of 2.55 for a given SFR; if 
the duration of star formation is decreased to 10 Myr 
(1 Myr), then the luminosity decreases by a factor of 
1.53 (12.9). 



4. APPLICATION TO (7-DROPOUT GALAXIES IN THE 

HDF 

4.1. Sample selection 

The above calibrations were applied to a sample 
of F300W (U) dropouts selected from the HDF V2.0 
catalog (Williams et al. 1996). We will compare the 
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Fig. 3. — The difference between photometric and spectroscopic plotted as a function of redshift. The x marks indicate the 
individual measurements of the 15 template spectra. The circles represent the median of the measurements and the errorbars the 
16th and 84th percentiles (i.e. m la errorbars). The solid line is the adopted polynomial calibration fit. 



results from this sample to two other U -dropout sam- 
ples extracted from the HDF by Madau and collab- 
orators - in M96 and Madau et al. (1998; hereafter 
M98). M96 lays out all selection criteria in detail and 
lists all galaxies selected, and in addition it is the pi- 
oneering work in this field. M98 improves on M96 
in several ways (see below). However, the results of 
M98 are model dependent, their [/-dropout selection 
criteria are not explicitly given, and their final sam- 
ple is not listed. Hence we adopt M96 as our primary 
comparison, and we model our sample selection on 
the criteria given therein. 

Our sample selection differs from that of M96, pri- 
marily because we use (^606 — ^814)ab rather than 
(-B450 — hi^AB as a color selection criterion. This 
is done because for {/-dropouts, the F450W filter 
is much more strongly affected by Lya forest ab- 
sorption than F606W. In our model, selecting by 
(^606 — ^814 )ab places a fairly strong limit on (3 and 
hence on ^i 6 oo- Figure | shows the (^06 - hu)AB, 
(^300 — -B45o)ab two color diagram, and our selec- 
tion box. Our selection criterion (V606 — ^814)ab < 
0.5 was chosen to encompass the strong [/-dropout 
plume but avoid relatively low-z interlopers. Our 
other selection criteria are virtually identical to M96; 
(^300 — -B45o)ab > 1.3 is the crucial [/-dropout se- 
lection criterion, while -B450 < 26.79 ABmag is the 
magnitude limit 7 . Only sources in the WF chips are 
selected. 



Another important difference with M96 is that 
their color selection box is not rectangular. They 
have an additional constraint ([/300 — -B45o)ab > 
(-6450- ^814)ab + 1-2, clipping the bottom right corner 
out of the box. This was introduced to avoid regions 
that model calculation indicate may be populated by 
low redshift galaxies. However, this corner, while rel- 
atively sparsely populated in the ([/300 — -B45o)ab ver- 
sus (-B450 — ^814) ab two color diagram, contains about 
Y 3 the galaxies in the HDF with measured 2 < z < 4 
(Dickinson 1998). Since the clipped corner is red in 
(-B450— hu)AB the clipping biases against large ^4i600- 
The color selection criteria of M98 are illustrated in 
Madau (1997), where they were revised in order to 
incorporate more galaxies with 2 < z < 4. Less of 
the bottom right corner is clipped than in M96. We 
have opted for a rectangular selection box in order 
to minimize this bias and keep the color selection as 
simple as possible, and yet be very close to the M96 
selection box. A total of 100 sources meet our selec- 
tion criteria. These are listed in Table |. Phot ome- 
try for these galaxies can be found in Williams et al. 
(1996). In comparison M96 selected 69 [/-dropouts, 
and Madau (1997) note that there are ~ 100 galaxies 
in their revised selection box. Dickinson (1998) using 
more robust photometry, a somewhat different color 
selection and a fainter limiting magnitude, isolates 
187 [/-dropouts in the HDF. 

A small difference with M96 concerns how sub- 



7 Nominally, we have also adopted Veoe < 28.0 ABmag in our selection. However, for our other selection criteria no objects 
approach this limit. 
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0.5 1 1.5 1 2 3 4 

(^606 ~~ ^BU^AB (^450 ~~ ^B^AB 

Fig. 4. — Two color diagrams of all objects (parents and daughters) in the HDF V2.0 catalog (Williams et al. 1996) with 
B450 < 26.79 ABmag. The dashed line in the left panel marks our adopted selection box for [/-dropouts, while that in the right 
panel marks the selection box of M96. The symbol correspondence is as follows: large triangles - spectroscopically confirmed 
galaxies with 2 < 2 < 4; squares - galaxies with spectroscopic z < 2; small triangles - other galaxies meeting our [/-dropout 
selection criteria (i.e. they must be in the box at left); dots - other objects not meeting our selection criteria. Some of these latter 
objects are nevertheless in our selection box. However they do not pass our simple parent-daughter separation algorithm, and are 
not included (and hence are not highlighted) in order to avoid double-counting. 



Table 3 

List of [/-dropouts selected in the HDF 



2-35.0 


2-689.2 


3-389.0 


4-317.0 


2-76.0 


2-726.0 


3-504.0 


4-363.0 


2-80.0 


2-741.0 


3-510.0 


4-382.0 


2-82.0 


2-751.0 


3-515.11 


4-389.0 


2-122.0 


2-793.0 


3-550.0 


4-445.0 


2-127.0 


2-802.11112 


3-593.0 


4-488.0 


2-239.0 


2-822.0 


3-633.1 


4-491.0 


2-321.2 


2-824.0 


3-674.0 


4-497.0 


2-373.0 


2-889.0 


3-677.0 


4-555.1 


2-381.12 


2-890.0 


3-736.0 


4-557.0 


2-392.0 


2-901.0 


3-748.0 


4-576.0 


2-449.0 


2-903.0 


3-813.0 


4-588.0 


2-454.0 


2-916.0 


3-875.0 


4-590.0 


2-456.12 


2-949.0 


3-902.0 


4-599.0 


2-456.22 


2-973.11 


3-916.0 


4-602.0 


2-525.0 


3-41.0 


4-52.0 


4-639.0 


2-547.0 


3-83.0 


4-83.0 


4-660.0 


2-565.0 


3-118.0 


4-85.2 


4-676.0 


2-585.1 


3-126.0 


4-96.0 


4-681.0 


2-591.0 


3-221.2 


4-109.0 


4-724.11 


2-594.0 


3-243.0 


4-194.0 


4-825.0 


2-604.1 


3-286.0 


4-252.0 


4-858.0 


2-637.0 


3-296.1 


4-274.0 


4-878.1 


2-643.0 


3-321.2 


4-289.0 


4-926.2 


2-689.11 


3-335.0 


4-316.0 


4-936.0 
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structure is included. The HDF catalog lists ob- 
jects found with the FOCAS software (Jarvis & 
Tyson 1981), which often splits a "parent" object into 
"daughters". Both parents and daughters are listed 
in HDF V2.0. M96 carefully evaluate the image of 
each candidate [/-dropout to determine whether it 
is likely to apply to the whole parent object or the 
daughter (s). Here, in cases where both daughter (s) 
and corresponding parent meet our selection crite- 
ria, only the parent object was included in our list of 
[/-dropouts. This could cause the integrated 1600A 
luminosity density pieoo to be overestimated. This is 
because not all daughters need satisfy the [/-dropout 
selection criteria for the parent to satisfy these same 
criteria. To estimate the maximum magnitude of 
this effect we considered the opposite extreme - when 
both parent and daughter(s) meet the selection crite- 
ria only the (lowest level) daughters that satisfy the 
selection criteria were included. This algorithm, re- 
sults in pi6oo being lower than our adopted sample 
by 16% (24% after absorption correction). This algo- 
rithm will underestimate P1600 because some daugh- 
ters, having the right z to be selected as [7-dropouts, 
will fail the magnitude limit, or because photomet- 
ric errors (larger for the fainter daughters) may shift 
some daughters out of the color selection box. The 
ideal parent-daughter separation will, therefore, lie 
somewhere between these two extremes, and hence 
the effect on pi6oo is likely to be a ~ 10% effect. 

4.2. Application of calibration 

The top panel of Figure ^ shows the color (/?) abso- 
lute magnitude M 1600,0 of our HDF U -dropout sam- 
ple. Here, (3 was calculated from the (V<306 — -^814)ab 
color listed in the HDF V2.0 catalog (Williams et al. 
1996) using eq. |l4|; where z was taken from the com- 
pilation of Dickinson (1998) where possible, other- 
wise we adopt an average value (z) = 2.75. This is 
the same [/-dropout (z) adopted by M96, and very 
close to {z) = 2.74 of the galaxies in our sample 
with measured z. The apparent AB magnitude at 
A = 1600A(1 + z) was interpolated between the V<306 
and the I814 magnitudes. This is effectively a k cor- 
rection. The y axis shows the absolute AB magnitude 
at Aq = 1600A corrected for UV absorption following 



eq. 11, while the broken diagonal lines indicate con- 
stant M1600 if there were no Aiqqq correction. The 
hatched region shows the region that is incomplete 
due to our selection criteria. The lower selection limit 
is slanted because of the ^4i600 correction and because 
of color terms in the k correction and in the transla- 
tion of the -B450 limiting magnitude into the Vqoq and 
/8i4 bands. 

The middle panel of Fig. |5| shows the (3 distribu- 
tion of the [/-dropouts, while the bottom panel shows 



that of the calibrating local sample for comparison. 
The color distribution of the local sample spans a 
larger range than the [/-dropouts. The two-sided 
Kolmogorov-Smirnov test gives a 0.3% chance that 
the two distributions were randomly drawn from the 
same parent distribution. These differences are par- 
tially due to our (V606 — /si4)ab < 0.5 selection limit 
which corresponds to /3 ^ —0.42 for z = 2.75, or 
equivalently ^4i600 ^ 3.6 mag. The local sample has 
no specific (3 selection. If we apply a (3 < —0.42 limit 
to the local sample then the two-sided Kolmogorov- 
Smirnov test gives a 3% chance that the local and U- 
dropout sample were randomly drawn from the same 
parent distribution; i.e. the two samples are more sim- 
ilar when similar (3 selection criteria are applied. The 
[/-dropout sample is bluer than the local sample hav- 
ing (3 < —0.42. The 20th percentile, median and 80th 
percentile of the [/-dropout distribution are -1.9, - 
1.6, and -1.2, while for the (f3 limited) local sam- 
ple they are -1.8, -1.4, and -0.8, respectively. One 
shouldn't read too much into this difference, in that 
the two samples were selected quite differently. The 
[/-dropouts have well defined magnitude and color 
limits in the rest-frame UV while the local sample is 
composed of galaxies observed with IUE for a wide 
range of programs with heterogeneous initial selection 
criteria. 

Luminosity functions at 1600 A, </>i6oo are shown 
for our U dropout sample in Fig. [6|. Two cases are 
shown: with (0i6oo,o) and without (0i6oo) the Aiqqq 
correction. The shape of </>i600 agrees well with that 
given by Dickinson (1998) for Mi 600 ^ -19 AB- 
mag. The drop-off in the last bin of our (fiieoo is due 
to incompleteness. Dickinson's luminosity function 
extends further because of his superior photometric 
technique and consequent fainter limiting magnitude. 
Our luminosity function is 0.3 dex higher than that 
displayed by Dickinson because he renormalizes to 
match the luminosity function of [/-dropouts selected 
from a larger ground based survey (Steidel et al. 
1998). This factor of two difference probably reflects 
the field to field variation in pi6oo- 

The 1600A luminosity functions extend well be- 
yond the range given for a sample of local UV se- 
lected galaxies by Buat & Burgarella (1998). Their 
luminosity function shows a distinct break at Lbol ~ 
10 10 L Q , corresponding to Mi6oo,o = —18.7 ABmag; 
a SFR of 4 M.q year -1 . The brightest galaxies in 
our [/-dropout sample have M1600 = —21.7 ABmag, 
Mieoofl = -24.1, corresponding to a SFR of 19, 220 
A^oyear -1 with and without absorption correction, 
respectively. This is further confirmation of the asser- 
tion that [/-dropouts are much more luminous than 
local UV selected starbursts. However, there are local 
starbursts with similarly high SFRs. These are the 
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Fig. 5. — The top panel shows the UV color j3 verus absorption corrected absolute magnitude Mi 60 o,o diagram for 
our sample of [/-dropouts. Large circles mark spectroscopically confirmed galaxies, while the rest of the [/-dropouts 
are marked with triangles. The broken lines mark equal M1500 without dust absorption correction. The hatched region 
marks the incomplete region according to our selection criteria. The slanted boundary is calculated for galaxies at the 
limiting -B450 = 26.79 ABmag with z — 2.75. It follows from eqs. |ll| and |l4] and relates the -B450 limiting magnitude 
to the Vqo6 and is 14 bands using (B450 — I&ia)ab = 0.19 + 2.09(V606 — ^814)ab, which was found using a linear least 
squares fit to the [/-dropouts having -B450 < 26 ABmag. The middle panel shows the (3 distribution of these galaxies. 
The solid histogram includes only the spectroscopically confirmed [/-dropouts. For comparison the f3 distribution of the 
local starburst sample (from Fig. |l|) is shown in the bottom panel. 
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Fig. 6. — Luminosity function (4>) at A = 1600A of our sample of [/-dropouts, with (top) and without dust absorption correc- 
tions. The error bars show the vJV uncertainties. The dotted line shows the completeness limit. Data to the left of these lines are 
complete, to the right they are incomplete. The top axis of each plot converts M1600 to star formation rate. The units for cj> are 
mag -1 Mpc -3 . The space density and luminosity of SCUBA sources, as estimated by Lilly et al. (1998), is indicated in the upper 
panel with a triangle. 
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ultra-luminous infrared galaxies (Sanders & Mirabel, 

1997) which are selected in the far infrared. 

From direct summation we derive the total ab- 
sorption corrected UV luminosity of our selected 
[/-dropouts to be ELi6oo,o = 3.8 x 10 34 ergs -1 . 
Assuming that the redshift range of z = 2.0 
to 3.5 is uniformly sampled, then the volume 
containing the [/-dropouts is 1.64 x 10 4 Mpc 3 
(M96), and the corresponding volume emissivity is 
,01600,0 = 1.4 ± 0.3 x 10 27 ergs" 1 Hz -1 Mpc" 3 ( PS fr = 
0.18 M® yr _1 Mpc -3 ). The error quoted here is the 
quadratic sum of the zeropoint errors from equations 
[10| , 11 and [14|, and the sampling error determined 
from the bootstrap technique (e.g. Babu & Feigel- 
son 1996). This pi600,0 is actually a lower limit for 
the volume sampled by the HDF because we have 
made no completeness correction. This raises a dim- 
cult issue. The standard practice is to estimate the 
completeness from the luminosity function, without 
dust absorption. Doing so, Madau (1997) estimates 
a completeness correction factor of 1.45. Figures ^ 
and H show that in the presence of dust, the limiting 
-^1600,0 is n °t uniform. Rather it is much shallower 
for highly reddened galaxies than for lightly reddened 
sources. Hence the completeness correction could be 
higher. In addition there may be field to field dif- 
ferences in /?i600,0) due to large scale clustering. Such 
variations could be as high as a factor of 2 (Dickinson, 

1998) . 

How accurate is our [/-dropouts selection criteria 
for selecting high-z galaxies? Our sample includes 
all but three of the spectroscopically confirmed HDF 
galaxies with z > 2 as compiled by Dickinson (1998). 
The galaxy 1-54.0 was missed because it appeared on 
the PCI chip rather than one of the WF chips; 3- 
512.0 has z = 4.02 and is -B-dropout and hence too 
red to be a [/-dropout; and 3-577.0 was not included 
because it is too faint for inclusion in our catalog. 
The remaining 23 spectroscopically confirmed high-z 
galaxies account for 47% of pi600,0 • We cross-checked 
our [/-dropout selection against a list of all redshift 
measurements in the HDF (kindly made available by 
Mark Dickinson). None of our selection corresponds 
to galaxies with known z < 2. Dickinson (1998) notes 
a 90% success rate in finding high-z galaxies in the 
appropriate redshift range using the U -dropout tech- 
nique on their ground based sample. While the pho- 
tometric bands and selection criteria are somewhat 
different from ours, this suggests that the contribu- 
tion of low-z interlopers to pi600,0 should not be much 
higher than 10%. 

4.3. The net effect of dust absorption 

In this paper we are primarily interested in the 
effects of dust on the high-z UV emissivity. If 



we calculate SL1600 and P1600 without correcting 
for dust absorption, the quantities are a factor 
of 5.4 ± 0.9 lower than the above estimates, e.g. 
Pieoo = 2.7 ± 0.2 x 10 26 ergs -1 Hz -1 Mpc -3 . This 
is the net effect of dust absorption on our sample 
of [/-dropouts. For comparison, M96 give pwoo = 
1.6 x 10 26 ergs -1 Hz -1 Mpc -3 for their sample of U- 
dropouts. The samples are identical in their limiting 
magnitudes, so the difference in the uncorrected emis- 
sivities is due to the sample selection. In particular, 
Fig. |4| shows that the portion of SL1600 missing from 
the M96 sample is contained within the corner clipped 
from their otherwise rectangular color selection box. 
This corner contains preferentially reddened galax- 
ies, hence this clipping biases against large Aiqqq. 
Relative to pimo estimated by M96 our emissivity 
is a factor of 9.0 larger. We reiterate that the actual 
absorption-corrected pi600,0 will be even higher than 
our estimates because we have not corrected for com- 
pleteness, and we are not sensitive to galaxies having 
-4i600 > 3.6 mag. 

More recently, Madau (1997) and M98 
now estimate an uncorrected pimo = 2.6 x 
10 26 ergs -1 Hz -1 Mpc -3 for the volume sampled by 
the HDF [/-dropouts, only 4% less than our esti- 
mate without dust correction. It is higher than in 
M96 because, as noted above, they adopt significant 
completeness corrections (a factor of 1.45) and use 
revised [/-dropout selection criteria that clip less of 
the bottom-right corner in the ([/300 — -B45o)ab vs. 
(£450 - /suUb plane. 

The 1600A dust absorption factor, 5.4±0.9 (equiv- 
alent to ^4i600 = 1-8 ± 0.2 mag), derived here for the 
[/-dropouts is lower than our original estimate of a 
factor of 15 (Paper II), albeit within the factor of 
three uncertainty we originally estimated. Sawicki 
& Yee (1997) using the Calzetti (1997) absorption 
law, derived similarly high dust absorption factors of 
15 to 20. Our new estimates are somewhat lower 
for a combination of reasons. Firstly, we now fit the 
IRX1600 - relation instead of adopting an absorp- 
tion law (Calzetti et al. 1994) that is unconstrained 
by FIR data. Secondly, in Paper II we adopted the 
median absorption of the spectroscopically confirmed 
[/-dropouts. Because of the Mi6oo,o — P correlation 
(Fig. ||) this biases the net absorption. Finally, when 
deriving from (Vgo6 — ^814)ab we now add a z de- 
pendent correction (eq. |l4|) to remove Lyman-forest 
and Lyman edge reddening, thus resulting in bluer 
colors and a smaller derived j4i600- 

The uncertainty in our new estimate of the mean 
dust absorption factor of [/-dropouts amounts to 
16%, and includes zeropoint and sampling errors. It 
is likely to underestimate the true uncertainty since 
it does not include field to field variations. M98 and 



15 



Pettini et al. (1998) estimate 1500A dust absorption 
factors of about three. Thus the difference between 
these groups and ours has narrowed considerabley, 
and is probably within realistic errors on the dust 
absorption factor. However we reiterate our caution: 
5.4 ± 0.9 represents a lower limit to the true mean 
1600 A dust absorption factor of galaxies at z ~ 3. 
As discussed in § |6.3| , true 1600A dust absorption fac- 
tors on the order of 10 are still not ruled out. 

Buat & Burgarella (1998) using a method very sim- 
ilar to ours, derive a somewhat lower mean absorption 
of 1.2 mag at A = 2000 A for a local sample of star- 
bursts selected to be visible in both the UV and FIR. 
The difference probably is largely due to the lower lu- 
minosity of the galaxies in their sample, which has a 
median L hoi ps 10 9 L e (SFR « 0.16 M Q year ). The 
[/-dropout sample has a median L^ Q \ = 8.5 x 10 10 L & 
(SFR 13A4© year -1 ), nearly a factor of 100 higher. 
Because of the Mi6oo,o — P correlation (Fig. ||) , we ex- 
pect the local sample to have a lower average amount 
of dust absorption (cf. Heckman et al. 1998). 

In summary, our value for the dust absorp- 
tion corrected UV emissivity preoo.o = 1.4 x 
10 27 ergs _1 Hz _1 Mpc~ 3 at (z) = 2.75 is a factor of 
9.0 higher than the uncorrected value given by M96 
and a factor of 5.4 times higher than the uncorrected 
value given by M98. The increased UV emissivity rel- 
ative to M96 is due to a combination of an improved 
[/-dropout selection (factor of 1.7), and correction for 
dust absorption (factor of 5.4). 

5. TESTS OF THE STARBURST HYPOTHESIS 

Our determination of the UV opacity in high-z 
galaxies relies on the similarity of the [/-dropouts to 
local starburst galaxies; the intrinsic SEDs of the two 
groups must be similar for our method to work. In 
§|2] we reviewed how the overall spectral properties of 
the two samples are similar. Here we concentrate on 
recent observational studies of high- z galaxies that 
are especially relevant to estimating UV opacities. 

5.1. Rest frame optical emission lines 

If the high- z galaxies are intrinsically redder in the 
UV than local starbursts (i.e. (3q > —2.2), then we 
would overestimate ^4i600 and hence pi600,0- This 
could occur if the high-z galaxies are dominated by 
older stellar populations (which would have a larger 
z of formation) or, if for some reason, their stellar 
IMF were deficient in the highest mass stars. Rest- 
frame optical emission lines present an ideal test for 
this hypothesis. Starbursts are by nature dominated 
by highly ionizing populations. Indeed the calibrat- 
ing IUE sample used in §|3| were classified as star- 
bursts primarily by the presence of intense narrow 
optical emission lines, arising in the ionized medium 



surrounding the starbursts. Hence, the ratio of line 
flux to flux density /1600 is a measure of the color 
from below the Lyman-limit to 1600A, modulated by 
dust absorption. 

For z = 2 — 3.5 rest-frame optical emission lines 
are redshifted into the observed-frame near infrared. 
With careful, and very long exposures, the bright- 
est optical lines are detectable in the brightest U- 
dropout galaxies with four meter class telescopes. 
Detections and fluxes have been reported by several 
groups including Pettini et al. (1998), Wright & Pet- 
tini (1998), Bechtold et al. (1997), and Bechtold et al. 
(1998). The detections correspond to [/-dropouts 
discovered from the ground or in the HDF, as well 
as other sources that meet the photometric criteria 
for being [/-dropouts. Figure [?] shows the line fluxes 
of these sources normalized by the rest frame 1600 
A flux density /1600 of the sources compared to (3. 
Both /1600 and were estimated from published pho- 
tometry, where possible (otherwise from published 
spectra) . The pseudo-equivalent widths inline) / /1600 
have been corrected for Galactic extinction and cos- 
mological stretching only. No attempt has been made 
to correct for internal absorption. Results for three 
different lines are shown - H/3, Ha, and [O ill] 5007A . 
In cases where the [O ill] 4959 A line was detected, it 
was multiplied by 3.0 to convert it to an [O m] 5007A 
flux. A sample of local starburst galaxies observed 
by Storchi-Bergmann et al. (1995) and McQuade, 
Calzetti, & Kinney (1995) is shown for comparison. 
The relevant data is listed in Table p]. This UV se- 
lected sample has a high degree of overlap with the 
sample used in S|| 

Figure shows that the detected [/-dropouts have 
line flux to UV flux density ratios consistent with 
local starburst galaxies. The one exception is the 
Ha observation of the lensed galaxy MS1512-cB58 by 
Bechtold et al. (1997). Wright k Pettini (1998) have 
also detected this galaxy in H/3 emission, and the cor- 
responding H/3 pseudo-equivalent width is well within 
the scatter of points for local starbursts as shown in 
the bottom panel of Fig |?]. The discrepancy between 
the Ha and H/3 measurements may be due to the 
lensed nature of the source, or perhaps may reflect the 
difficulty of the Ha observations which were obtained 
at A = 2.443 /im, the edge of the K band window. At 
this stage, none of the groups report non-detections 
in emission lines of known [/-dropouts. We conclude 
that [/-dropouts are ionizing sources to the same de- 
gree as local starbursts, and hence do not have an 
intrinsically redder continuum than the local sample. 

5.2. Far infrared and radio continuum observations 

If there is significant UV extinction of galaxies at 
z ~ 3, then perhaps the dust can be detected directly 
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Fig. 7. — Line flux to UV flux density ratio compared to UV spectral slope. The open circles show a local sample of UV 
selected starbursts, using data originally published by Storchi-Bergman et al. (1995) and McQuade et al. (1995). The labelled 
triangles represent seven (7-dropouts observed by Pettini et al. (1998), Wright & Pettini (1998); Bechtold et al. (1997) and Bechtold 
et al.(1998). The label - galaxy name correspondence is (1) MS1512-cB58; (2) Q0000-D6; (3) Q0201-C6; (4) Q0201-B13; (5) B2 
0902-C6; (6) DSF 2237a-C2; and (7) A2218-A384. These data have been corrected for the redshift stretching of the spectral energy 
distribution. 



in the rest frame FIR or indirectly in the radio via the 
radio-FIR correlation (Helou 1991). Two groups have 
recently claimed (marginal) detections of sources in 
and near the HDF in these wavelength regimes some 
of which which they attribute to z ~ 3 star forming 
galaxies. Here we show that the number and fluxes 
of these marginal detections are consistent with the 
expectations of our dust absorption model. 

The sub-mm regime, corresponding to the rest 
frame FIR, is a promising window to look for high-z 
galaxies since the "negative ^-corrections" for dust 
emission can cancel out the effects of fade-out with 
distance, making dusty galaxies easier to detect 
with increasing z (Guiderdoni et al. 1997). Several 
groups have now obtained deep sub-mm imaging with 
SCUBA 8 on the 15m JCMT 9 (e.g. Smail, Ivison & 
Blain, 1997; Barger et al. 1998; Eales et al. 1998; 
Hughes et al. 1998; Lilly et al. 1998). All report 
850/xm detections at the mJy level that they largely 
attribute to galaxies having z > 1. Indeed one of 
the SCUBA detections, SMM02399-0136, has a clear 
optical counter-part and spectroscopic determination 
of z = 2.8 (Ivison et al. 1998). However, Lilly et al. 
(1998) note that a sizable fraction of these sources 
may have z < 1. Furthermore, the emission line spec- 
trum of SMM02399-0136 is indicative of a Seyfert 2 
galaxy, cautioning us that not all of the dust seen by 

8 Sub-mm Common User Bolometer Array 

9 James Clerk Maxwell Telescope 



SCUBA need be warmed by star- formation alone. 

The study of Hughes et al. (1998) is particularly 
relevant because it is of a field centered on the HDF. 
They report five sources above their detection limit 
of 2mJy. Since this limit is set by source confusion 
(beam FWHM~ 15"), none of their optical identifica- 
tions is totally secure. They attempt to improve the 
accuracy of their identifications using photometry, in- 
cluding photometric redshifts. They claim that 4/5 of 
their detections are best identified with sources hav- 
ing z ^ 2. The brightest of these has fsbo^m = 7 ± 2 
mJy. 

In Table || we present /sso^m estimates for HDF 
?7-dropouts based on our dust absorption model, and 
the optical photometry of Williams et al. (1996). To 
do this we used eqns. [K] and [Ll] to estimate the rest 
frame -Ffir from the rest frame -F1600 an d (3 estimates. 
-FpiR was then converted to /850/^m using the FIR SED 
of the prototypical IR luminous starburst Arp 220 as 
modeled by a spline fit through the the measurements 
of Rigopoulou et al. (1996) and Klaas et al. (1997). 
Only the galaxies with the brightest predicted /s50/xm 
values are listed. The table is arranged in order of 
decreasing predicted /sso^m- The corresponding rest 
frame 1600A flux densities (uncorrected for dust ab- 
sorption) do not come out to be similarly sorted. This 
is because it is the dust reprocessed flux that is im- 
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portant. Since our adopted [/-dropout sample is color 
limited (i.e. Aiqqq <; 3.6 mag), the brightest predicted 
infrared fluxes correspond to the brightest and red- 
dest {/-dropouts. The UV absorption ^4i6oo listed in 
this table is calculated from the (Vqoq — Jsi4)ab as 
outlined in 



3.3 



None of these galaxies have predicted /sso^m > 2 
mJy, consistent with the fact that none were detected 
by Hughes et al.. Likewise, none of the claimed opti- 
cal counterparts of /sso^m detections meet our criteria 
of being [/-dropouts. One is clearly a z ~ 1 galaxy. 
Three are fainter than our adopted inclusion limit. 
The remaining galaxy HDF850.3 is in the PCI chip 
and identified as 1-34.2. Its (Vqo6 — -^814 )ab = 0.33 is 
within our selection limits, but its ([/300 — -B45o)ab = 
1.08 color is marginally too blue to be considered aU- 
dropout, but consistent with a photometric z ~ 1.95 
estimated by Mobasher et al. (1996). 

The non-intersection of our [/-dropout catalogue 
and the 850/um sample of Hughes et al. suggests that 
the sub-mm detections correspond to infrared lumi- 
nous starbursts or AGN that are so dust enshrouded 
that they are not selected as [/-dropouts either be- 
cause they are too extincted and/or because they 
are too reddened. The Hughes et al. identification 
of two 850/xm sources (HDF850.1, HDF850.2) with 
fairly blue galaxies ((V<306 — -[si4)ab < 0.5) also sug- 
gests that the /3-IRXi6oo calibration breaks down for 
heavily extincted infrared-luminous starbursts. Fur- 
ther work is needed to confirm these identifications 
and secure reliable redshifts. 

The remarkably tight radio - FIR correlation 
(Helou 1991) provides an indirect means to probe 
the FIR emission, and hence dust emission of star- 
bursts. It is best defined at v = 1.4 GHz (21cm; 
e.g. Sanders & Mirabel, 1997). The radio emission 
of a wide range of star forming galaxies, including 
starbursts, is dominated by synchrotron emission at 
this frequency and up to ~ 30 GHz. Hence the cor- 
relation can be used to infer FIR fluxes from radio 
observations down to rest wavelengths of about 1 
cm. Do we expect this correlation to hold at high- 21 ? 
Inverse Compton scattering from cosmic microwave 
background photons provides an alternative mecha- 
nism to synchrotron emission for cooling the relativis- 
tic electrons. The energy density of these photons 
will scale as (1 + z) 4 . However it is the total en- 
ergy density of photons that is important for inverse 
Compton cooling. Even at z ~ 3 the energy density 
of cosmic microwave background photons will be an 
order of magnitude lower than internally produced 
energy density of starbursts having typical bolomet- 
ric surface brightnesses of S e = 3 x 10 10 L kpc~ 2 
(e.g. Paper II) at z = 3. Since synchrotron emission 
clearly occurs in the face of intense internally pro- 



duced emission in low-z galaxies, it does not seem 
likely that inverse Compton cooling will quench it for 
high- 2: galaxies where the total energy density is only 
marginally higher. 

Richards et al. (1998) present a deep imaging sur- 
vey of the HDF and surrounding areas obtained with 
the VLA at 3.5 cm (v = 8.46 GHz). Their 5a detec- 
tion limit is 9 /xJy, and their 3.5cr marginal detection 
limit corresponds to 6.3 /iJy. They marginally detect 
two HDF [/-dropouts, 4-52.0 and 2-449.0 as listed in 
the HDF V2.0 catalog. While the detections are pro- 
visional, Richards et al. are confident that they are 
not spurious because they also are detected in a 1.4 
GHz image of the HDF. 

The logarithmic ratio of f u flux densities in the rest 
frame radio (1cm = 30 GHz) and UV (1600A = 1.9 
PHz) and can be written as: 



log 
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where q v is the radio-FIR logarithmic ratio: 
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At A = 21 cm, q v {2lcm) = 2.35 (Sanders & Mirabel, 
1997). Since starbursts typically have radio spec- 
tral slopes a » —0.7 (where f u oc u a ), we ex- 
pect q v (\cm) = 3.26. For the prototypical starburst 
galaxy M82, gj,(lcm) = 3.16 is derived from the data 
listed by Carlstrom & Kronberg (1991). This is some- 
what lower than expectations because free-free emis- 
sion is starting to make a contribution to the lcm flux 
of M82 (Carlstrom & Kronberg 1991). We adopt a 
rest frame ^(lcm) = 3.2 and calculate the expected 
f u at 3.5cm (/3.5cm) of the [/-dropouts in our HDF 
sample using the above formulation, and modest (fac- 
tor of a few at most) u k n corrections to transfer from 
the rest frame wavelengths of 1600A and lcm, and 
the observed frame /si4 an d A = 3.5 cm observations. 
Table || lists the results of these calculations for the 
expected 10 sub-mm brightest [/-dropouts. Note that 
the f u ranking of the radio predictions is different 
from that of the FIR predictions. This is because the 
k corrections in the two wavelength regimes work in 
opposite direction with redshift at z ~ 3. 

The sources in Table [| include the two detections of 
Richards et al. (1998). Note that all of the galaxies in 
this table have predicted /3.5cm below the formal de- 
tection limits of Richards et al. . However the top five 
/3.5cm estimates are not far below their marginal de- 
tection limit, and include these two sources. In such 
cases we expect the fortuitous coincidences of uncer- 
tainties to allow some positive detections. These un- 
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Table 4 

hdf [/-dropout fir and radio flux predictions 











A: 


(z + l)1600A 


850Atm 


3.5cm 


3.5cm 












Observed 


Predicted 


Predicted 


Observed 


Name 


2 


^814 


P 


^1600 


fu 


fu 


fu 


fu 






[ABmag] 




[mag] 


My] 


[mJy] 


My] 


My] 


4-52. a 


2.931 


24.06 


-0.66 


3.05 


0.63 


1.8 


3.7 


6.3 ±2.6 


4-878.1 


(2.75) 


23.28 


-1.06 


2.32 


1.35 


1.6 


3.7 




4-639.0 


2.591 


24.16 


-0.73 


2.92 


0.52 


1.1 


2.8 




2-585.1 


2.002 


23.50 


-0.73 


2.93 


0.76 


1.1 


4.7 




4-555.1 


2.803 


23.33 


-1.27 


1.93 


1.39 


1.1 


2.3 




4-445.0 


2.268 


23.75 


-0.80 


2.79 


0.69 


1.0 


3.5 




4-382.0 


(2.75) 


24.10 


-0.95 


2.52 


0.61 


0.9 


2.2 




2-127.0 


(2.75) 


24.72 


-0.79 


2.81 


0.33 


0.7 


1.6 




2-449.0 b 


2.008 


23.42 


-1.14 


2.17 


1.01 


0.6 


2.7 


7.0 ±2.0 


4-858.0 


3.216 


24.06 


-1.40 


1.70 


0.78 


0.6 


1.0 





a Source 3647+1255 in Richards et al. (1998) 
b Source 3648±1416 in Richards et al. (1998) 



certainties include the measurement bias mentioned 
above, the flux uncertainties listed by Richards et al. 
(~0.1 dex uncertainty), the intrinsic scatter about the 
radio-FIR relationship (~ 0.3 dex uncertainty; Dev- 
ereaux & Eales 1989), and the random scatter about 
the ^1600 versus /3 fit (~0.1 dex uncertainty). The to- 
tal random uncertainty is ~ 0.3 dex or a factor of ~ 2. 
Hence random errors alone can account for the differ- 
ence in the predicted and observed /3.5cm in the two 
galaxies detected by Richards et al.. Why aren't these 
galaxies also detected in the rest frame FIR? Scaling 
our predictions by the /3.5cm observed by Richards 
et al., we then expect /sso^m = 3.2 ± 1.3, 1.6 ± 0.5 
mJy for 4-52.0 and 2-449.0 respectively. The latter is 
below the Hughes et al. detection limit, while the for- 
mer is above it but by less than the flux uncertainty, 
and so could easily be missed. We conclude that the 
Richards et al. radio detections are consistent with 
the Hughes et al. SCUBA non-detections. 

Our calculations show that [/-dropouts in the HDF 
are just at the detection limit of state of the art deep 
rest frame FIR and radio observations. Somewhat 
deeper observations in these bands of either the HDF 
or other fields containing [/-dropouts (preferably the 
brightest and reddest ones) have the potential for di- 
rectly testing the - IRX1600 relationship and ver- 
ifying the existence of the radio-FIR correlation at 
z rj 3. Hence they could be used to directly test 
the form of the UV absorption law. For example, 
the SMC law absorption model adopted by Pettini 
et al. (1997) would predict /sso^m < 0.4 mJy, and 
/3.5cm < l^Jy f° r an the [/-dropouts in our sample. 
Hence mJy level imaging at 850/um and /uJy level ra- 
dio continuum imaging should detect the brightest 



[/-dropouts if our absorption model is correct, but 
would not detect any if lightly absorbed models were 
correct. 

6. DISCUSSION 

6.1. Color-luminosity correlation 

The color-magnitude diagram of HDF [/-dropouts, 
Fig. |5], shows that the highest luminosity galaxies, 
with or without A1600 correction, tend to be red. 
Similarly the bluest galaxies tend to have the low- 
est luminosity. That is, there is a color-luminosity 
correlation (also noted by Dickinson, 1998, private 
communications). The (3 - luminosity correlation is 
similar to that seen in local starburst galaxies (Heck- 
man et al. 1998). As is the case for local starbursts, 
this correlation is not solely due to selection effects. 
While Fig. || clearly shows that we select against 
low-luminosity red [/-dropouts, there is no selection 
against high luminosity very blue objects. The more 
luminous starbursts have a higher dust content. Pre- 
sumably, dust content is related to metallicity. This 
is certainly born out for local starbursts where [O/H] 
correlates with (3 (Heckman et al. 1998). From Paper 
II, we know that the [/-dropouts have similar surface 
brightnesses, hence the most luminous ones are the 
largest, and presumably the most massive galaxies. 
Putting these together implies a mass-metallicity re- 
lationship for [/-dropouts like that seen in starburst 
(Heckman et al. 1998) and normal galaxies in the lo- 
cal universe. The presence of dust in the most lumi- 
nous [/-dropouts implies that they were previously 
enriched with metals. Hence we are not seeing the 
most luminous galaxies at the instant they "turn on" . 
Their first generation of star formation must have oc- 
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curred at z > 2.75. 

6.2. Optical line to UV continuum ratios as 
reddening indicators 

Figure [7| illustrates the (un) usefulness of line to 
UV continuum ratios as indicators of dust absorp- 
tion; the y axis is proportional to the ratio of the 
SFRs one would derive from emission lines and the 
UV continuum respectively, in the absence of dust 
absorption. Pettini et al. (1998) argue that the sim- 
ilarity of the SFR derived from H/3 fluxes and the 
UV continuum indicates that dust absorption is not 
extreme. However, Fig. shows little correlation be- 
tween Fnp/fiQQQ with (3. In local starburst systems, 
we know that (5 is a good dust absorption indicator 
(Fig. |); apparently F Rf3 /f 1600 is not. The only ra- 
tio that noticeably correlates (positively) with f3 is 
-^H«//i600- However, the slope of the correlation is 
about half that expected by application of the LMC 
extinction law. The Fro m] //i600 ratio even shows 
an anti-correlation with (3. This is probably due to 
the correlation of [O/H] with (5 found for starbursts 
(Heckman et al., 1998) - highly reddened starbursts 
have higher metallicities, lower electron temperatures 
in the ionized medium, and hence lower [Oiii]/H/3 ra- 
tios than low reddening starbursts. 

The .FW/ieoo and F H p/f imo results illustrate 
a result that is often overlooked: the dust ab- 
sorption estimated from continuum observations is 
usually less than that estimated from Hll emis- 
sion. This was first noted by Fanelli, O'Connell, & 
Thuan (1988) and subsequently confirmed by oth- 
ers including Mas-Hesse, Arnault, & Kunth (1989), 
and Calzetti et al. (1994). The lack of correla- 
tion between -FH/3//1600 an d indicates that the 
effective extinction law of starbursts is such that 
^4i600 (continuum) ss Ah/^Hii). While uncorrected 
UV continuum fluxes and H/3 line fluxes may imply 
the same SFRs, both are affected by similar amounts 
of dust absorption. Hence, emission line flux to UV 
continuum ratios do not provide strong constraints 
on dust absorption. 

6.3. SCUBA detections and the nature of 
U -dropouts 

The SFR of the various SCUBA detections is typ- 
ically derived to be > lOCMoyr -1 (e.g. Lilly et al. 
1998). As shown in Fig. ^, this is consistent with 
the SFR derived for the brightest and reddest U- 
dropouts, after dust absorption correction. Such high 
values for the SFR are unprecedented for local UV 
selected starbursts, and typically only seen in FIR 
selected starbursts such as Arp220. This suggests 
that the most luminous starbursts at z ~ 3 are less 
dust enshrouded than those seen locally - the ultra- 



luminous infrared galaxies. However, this conjecture 
is difficult to test, since there are very few UV obser- 
vations of ultra-luminous infrared galaxies (cf. Tren- 
tham et al. 1999), hence the amount of UV absorption 
they suffer has not been properly characterized. Also 
shown in Fig. || is Lilly's et al. (1998) estimate of the 
number density of the SCUBA sources in the red- 
shift range 2 < z < 3. This is also in fair agreement 
with the number density of the [/-dropouts having 
the same luminosity. This comparison suggests that 
the SCUBA sources may be largely [/-dropouts with 
modest amounts of dust absorption, rather than be- 
ing completely opaque in the UV. Furthermore, it 
suggests that much of the SFR density at z « 3 can be 
recovered from rest-frame UV measurements alone. 

Working against this suggestion is the fact that 
none of the HDF SCUBA sources detected by Hughes 
et al. (1998) correspond to known [/-dropouts. How- 
ever, the numbers are small, and the association with 
optical galaxies is still not secure. Hughes et al. es- 
timate psfr ~ 0.11 M.q yr _1 Mpc -3 at z 3 from 
their SCUBA detections which is about half of the 
Psfr we derive. This is only ~ five times larger than 
the M96 estimate (~ 3 times higher than the M98 
estimate). Like us they do not extrapolate the lumi- 
nosity function of their detections, hence the ^sfr is 
a lower limit consistent with our estimates. Alterna- 
tively, if the SCUBA sources are at z ~ 3 and a non- 
intersecting set with respect to [/-dropouts (as sug- 
gested by the Hughes et al. identifications), then the 
sum of the Hughes et al. /?sfr and that derived from 
our [/-dropout sample provides an improved lower 
limit: psfr ^> 0.30 M. & yr _1 Mpc~ 3 at z 3 which is 
a factor ^15 larger than that estimated by M96, and 
;> 9 times the uncorrected value estimated by M98. 

M98 argue against such a high value of /?sfr since 
it would over-produce low mass stars in the present 
universe if it went on for more than about one Gyr. 
However, this argument assumes that the Lyman- 
dropouts have an IMF with constant Salpeter (1955) 
slope over the mass range of 0.1 to 125 Mq. Re- 
cent observations of local starbursts suggest that their 
IMF is deficient in stars less massive than 3Mq with 
respect to the Salpeter IMF (Nota et al. 1998; Siriani 
et al. 1999) 

Since psfr is very model dependent, a more ro- 
bust prediction of interest is the lower limit to 
the FIR emissivity. Summing the predicted FIR 
fluxes of all the HDF [/-dropouts (using eq. 10) 
and normalizing by the volume yields pfir > 7.2 x 
10 29 ergs -1 Hz _1 Mpc~ 3 at z = 2.75 evaluated at 
a rest wavelength of 80/iin. If the SCUBA detec- 
tions are too reddened or dimmed to be selected 
as [/-dropouts then we can add them to our sam- 
ple to get an improved lower limit: pfir > 12 x 
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10 29 erg s _1 Hz^Mpc -3 . These are still lower lim- 
its because completeness corrections have not been 
made. These predictions can be tested directly with 
improved FIR, sub-mm, and radio observations (e.g. 
with the upcoming Sub-mm Array, and the SIRTF, 
FIRST, and PLANCK space missions). 

6.4. Application to other redshifts 

Our calibration of UV extinction is readily adapted 
to observations at redshifts other than z ~ 3 and 
hence other spectral domains. The main requirement 
is that the data include fluxes in two broad bands or 
coarse spectra covering the rest frame UV. Redshift 
information is also helpful, and that can be obtained 
photometrically to adequate accuracy from comple- 
mentary rest frame optical - near IR data (e.g. from 
the Sloan Digital Sky Survey, or deep pointed ob- 
servations). For example, serendiptous detections of 
galaxies with z ;> 5 (e.g. Franx et al. 1997; Dey et al. 
1998) indicate that galaxy formation had commenced 
by this epoch. One of these galaxies is a lensed arc 
at z = 4.92, whose broad band (observed frame) op- 
tical - near IR photometry yield (3 = —1.63 (Soifer 
et al. 1998), implying Ai6oo ~ 1-2 mag. Upcoming 
all sky surveys in the UV by GALEX (Martin et al. 
1997) and TAUVEX (Brosch, et al. 1994) will probe 
the cosmological UV emissivity to z <^ 2. Extending 
deep imaging with the HST into the near IR with 
NICMOS allows Lyman edge systems to be detected 
out to z ~ 8 (/814-dropouts). 

There are a few caveats to using our methods, espe- 
cially at \ow-z. Firstly, it is based on the assumption 
that the UV emissivity is dominated by starbursts. 
Possible contaminants are normal disk galaxies and 
old stellar populations. Ideally we want to include 
the UV light from normal disk galaxies when estimat- 
ing star SFRs. The problem is that it is not known 
whether normal disks follow the starburst reddening 
law derived here. Their cool and warm ISM is dom- 
inated by a disk morphology, while starbursts typi- 
cally have significant quantities of extra-planar ISM 
ejected by galactic winds (Heckman, Armus, & Miley 
1990). In Paper I we suggested that winds shape the 
ISM around starbursts into a foreground screen, thus 
producing the IRX1600 - P correlation. Further work 
has to be done to determine if normal disks have a 
reddening - UV absorption relation and to calibrate 
it. Fortunately normal disk and starburst galaxies 
can be distinguished in the UV by effective surface 
brightness (Paper II). Cosmological surface bright- 
ness dimming renders normal spiral galaxies difficult 
to detect by z 3 (Giavalisco et al. 1996; Hibbard & 
Vacca, 1997), hence they do not contaminate the U- 
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dropout sample. Old stellar populations produce sig- 
nificant UV "upturn" emission. If they significantly 
contaminate a UV survey then the assumption that 
UV light is due to recent star formation is seriously 
compromised. Old stellar populations can be distin- 
guished and removed from UV surveys by extending 
the broad-band SED into the rest frame optical. 

Secondly, our technique has been calibrated to a 
fairly modest A1600 ~ 5 mag (Fig. [|) . While the local 
calibrating sample contains several well known dusty 
starbursts (e.g. NGC 1614, NGC 3504, NGC 3256, 
NGC 7552), it doesn't contain any of the "Ultra- 
Luminous Infrared Galaxies" (ULIRGs) which are 
the most luminous starbursts known (Sanders & 
Mirabel, 1997). It is not known if they fall on or con- 
tinue the IRX1600 - P relationship since the number 
of these galaxies observed in the UV is rather mea- 
ger (Trentham et al. 1999). The amount of totally 
obscured star formation in the local universe can be 
estimated by applying our absorption algorithm to a 
blind UV survey such as will be produced by GALEX, 
after correction for contaminants as suggested above. 
This will yield an estimate of /0i6oo,o an d also the cor- 
responding density of UV flux processed to the FIR. 
Comparison with pfir derived from IRAS, ISO 10 , or 
SIRTF 11 observations will yield the fraction of com- 
pletely buried star formation. 

7. SUMMARY 

The strong similarity between local starburst galax- 
ies and Lyman limit systems at z ~ 3 provides a 
strong impetus to use starbursts as templates for in- 
terpreting the properties of high- z galaxies. The fact 
that starburst galaxies redden as more of their UV 
light is absorbed by dust is demonstrated by a strong 
correlation between their FIR to UV flux ratios and 
their ultraviolet spectral slopes j3 (color). Our A1600 
versus (3 calibration is a simple empirical fit to this 
correlation. Hence, by design it recovers the UV light 
that has been reprocessed to the FIR. 

We applied this calibration to a [/-dropout sam- 
ple selected from the HDF V2.0 catalog. This sam- 
ple has an identical limiting mag to that of Madau 
et al. (1996: M96). However, the color selection 
is somewhat different, resulting in a fairly uniform 
limiting A1600 ^ 3.4 mag. In contrast the M96 
[/-dropout sample has a non-uniform bias against 
large ^4i600- For our adopted cosmological param- 
eters of Hq = 50 kms -1 Mpc -1 , qo = 0.5, the 
rest-frame comoving UV emissivity of our sample 
is P1600 = 2.7 x 10 26 ergcm~ 2 s _1 Hz -1 Mpc -3 with- 
out applying ^4i6oo corrections. This is 1.7 times 
higher than the value derived by M96, which illus- 
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trates the effects of the differing color selections. Af- 
ter correcting for dust absorption we find pwoofl = 
1.4 x 10 27 erg cm~ 2 s _1 Mpc -3 , which is factors of nine 
and five times higher than the uncorrected values of 
M96 and M98, respectively. However, this pi600,0 is 
still a lower limit because we have not corrected for 
incompleteness, and we are not sensitive to galaxies 
having Aiqoq > 3.6 mag. 

This study and other recent observations further 
demonstrate the close correspondence between local 
starbursts and high-z Lyman limit systems. Firstly, 
we find a correlation between (3 (color) and Mi6oo 
(UV luminosity) for the [/-dropouts, similar to the 
correlation seen in local starburst galaxies (Heckman 
et al. 1998). To first order the UV luminosity gives 
the mass of recently formed stars, while /3 correlates 
well with metal content in local galaxies (Heckman 
et al. 1998). This implies that there is a mass metal- 
licity relationship out to z ~ 3. Since the highest 
luminosity galaxies are clearly reddened they must 
have produced enough metals to create the dust prior 
to the epoch of observation. Hence the first epoch 
of star formation occurred at z > 2.75 in the most 
luminous [/-dropouts. Secondly, observation of emis- 
sion lines in the rest frame optical (observed in the 
near-IR) show that [/-dropouts contain ionizing stel- 
lar populations to the same degree as local starbursts. 
This can be seen by comparing the line flux to UV 
flux density ratios with (5 for both samples. How- 
ever, neither the local nor the [/-dropout samples 
show a strong correlations in these diagrams as would 
be expected from naive application of UV extinc- 
tion/absorption laws. This suggests that the column 
density of dust towards Hll emission is higher than 
towards the UV continuum both in local starbursts 
and out to z 3. Consequently, the line flux to UV 
flux density ratio is not a good indicator of dust ab- 



sorption. Thirdly, the marginal detection of two HDF 
[/-dropouts at radio rest Ao ~ 1 cm by Richards et al. 
(1998) and non-detections of the [/-dropouts in the 
HDF by Hughes et al. (1998) is consistent with es- 
timates of the amount of UV flux reprocessed into 
the FIR. Deeper observations in both domains (by a 
factor of ~ \/2)) are required to test the dust repro- 
cessing model proposed here, and to confirm whether 
the FIR-radio correlation holds at z ~ 3. 

The methodology presented here is readily adapted 
to other redshifts. For those sources affected by mod- 
est amounts of absorption (^4i6oo ^ 5 mag), the in- 
trinsic UV flux can be recovered from UV quantities 
alone: flux measurements in two rest frame vacuum- 
UV bands are sufficient. These observations can also 
be used to estimate the total UV flux reprocessed 
into the FIR. Hence, lower limits to the intrinsic 
pxjVfi(z) and pfir(z) can be determined. Integration 
of Pfir(-z) and comparison to the FIR background 
spectrum (Fixsen et al. 1998; Hauser et al. 1998) will 
constrain the amount of star formation totally ob- 
scured by dust from the ultraviolet. 
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